Propagating phase boundaries as sonic horizons 
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If certain conditions are met, a propagating phase boundary can be a sonic horizon. Sonic Hawking 
radiation from such a phase boundary is expected in the quantum theory. The Hawking temperature 
for typical values of system parameters can be as large as ~ 0.04 K. Since the setup does not require 
the physical transport of material, it evades the seemingly insurmountable difficulties of the usual 
proposals to create a sonic horizon in which fluid is required to flow at supersonic speeds. Issues that 
are likely to present difficulties that are particular to this setup are discussed. Hawking evaporation 
of the sonic horizon is also expected and is predicted to lead to a deceleration of the phase boundary. 
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Condensed matter analogs of gravitational systems 
have brought purely theoretical ideas on very large-scale 
gravitational structures to almost within the realm of ex- 
perimental tests. A key breakthrough was the realization 
by Unruh j3| that black hole event horizons can be sim- 
ulated by sonic event horizons around "dumbholes" in 
fluid flow. Furthermore, theoretical results from the ap- 
plication of quantum fleld theory in curved spacetimes, 
such as Hawking radiation Q, should also apply to con- 
densed matter analogs, thus opening the door to novel 
experimental tests of gravitational phenomena. 

The fluid analog of black holes occurs when sub-sonic 
fluid flow upstream changes to supersonic flow down- 
stream. Now the fluid flow downstream is too rapid to 
allow sonic perturbations to flow upstream. Sound can 
travel downstream from upstream but not upstream from 
downstream, resulting in a "sonic horizon" at the loca- 
tion where the flow velocity equals the sound velocity. To 
the upstream observer, the downstream region is a hole 
from which no sound can emanate, that is, it is a dumb- 
hole. This is the classical picture. As in the black hole 
case, "dumbholes ain't so dumb" and quantization of the 
acoustic perturbations leads to sonic Hawking radiation 
from downstream to upstream. Unruh found that the 
sonic Hawking temperature is given by: 
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where fi is Planck's constant divided by 27r, ks is Boltz- 
man's constant, r is the radial coordinate, v = |v|, and 
the flow velocity v is assumed to be radially inward. The 
derivative of the velocity is evaluated at the location of 
the sonic horizon determined hy v — Cs where Cs is the 
sound velocity. A numerical estimate gives: 
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where the velocity gradient is taken to be Cg/R- The 
estimate clearly demonstrates the challenge: even when 



the fluid undergoes tremendous acceleration - velocity 
change from to 300 m/s within 1 mm or, equivalently, 
10^5 where g is the terrestrial acceleration due to gravity 
- the sonic Hawking temperature is extremely small. 

The immense practical problems in creating a sonic 
black hole and observing Hawking radiation have been 
described by Unruh leading to a very grim picture. 
The crux of the problem is that the flow of fluids and su- 
perfluids suffers from instabilities at the tremendous ac- 
celerations that are required for creating a sonic horizon 
with a measurable Hawking temperature. If one chooses 
to work with a system with low sound velocity, the insta- 
bilities can be evaded, but then the Hawking tempera- 
ture is unmeasurably low. At the moment these difficul- 
ties seem to be insurmountable, except possibly in Bose- 
Einstein condensates 0, |^ . (Proposals to test quantum 
processes in the presence of cosmological horizons in BEC 
analogs also appear promising ^6,].) 

For the idea to be described in the present paper, a 
generalization of Eq. made by Visser will be very 
important. Visser considered the case when the sound 
velocity itself is a function of space and time and ob- 
tained: 
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This seemingly small step opens up a conceptually new 
set of possible dumbhole realizations because now it is not 
necessary to accelerate the fluid; instead one can manip- 
ulate the function Cs(i,x). 

Along these lines, the proposal of Jacobson and Volovik 
(Sj (also see :9;]) employs a domain wall in '^He-A to pro- 
duce a variation of the "sound" velocity on micro- 
scopic scales (~ 500 A). Modest velocities of the domain 
wall lead to an enhanced Hawking temperature. The gain 
is however somewhat offset by the rather low sound ve- 
locity 3 cm/s). The Hawking temperature in their 
setup is 5/xK, whereas reliable thermometry in '^He cur- 
rently only goes down to 100/xK. Furthermore, gradients 
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FIG. 1: Essentials of the experimental setup described in the 
text. The system is kept at constant temperature and pres- 
sure. 



in the order parameter associated with the presence of 
the domain wall lead to additional sources of radiation 
that are not related to the sonic horizon. It remains to be 
seen if the Hawking radiation can be disentangled from 
these other effects. 

We envisage a setup somewhat related to the Jacobson- 
Volovik proposal and pictured in Fig. ^ The container 
is initially filled with a metastable phase (Phase 1) of a 
material. Suppose that Phase 2 is the stable phase under 
the conditions of the experiment. Then a perturbation 
will trigger bubble nucleation of Phase 2 and this bubble 
will grow. The velocity v with which the bubble grows 
will be subsonic in Phase 1: v < Ci where ci is the sound 
velocity in Phase 1 . The sound velocity C2 in Phase 2 will 
be different from that in Phase 1. We assume C2 < Ci. 
More crucially we take C2 < v. For the purposes of this 
paper, we will simply assume that the conditions (ci > 
V > C2) are met in the system we have chosen. Then 
sound emanating from within Phase 2 cannot catch up 
with the phase boundary. Thus the phase boundary will 
be a sonic horizon. By the usual arguments J7|, sonic 
Hawking radiation should be seen to emanate from the 
bubble [nl. 

The advantage of this setup is that it completely elim- 
inates any instabilities associated with supersonic flow. 
In certain situations, the bubble wall may push on the 
material outside. If the system is kept at constant pres- 
sure, this would cause some expansion of the occupied 
volume. The height occupied by the material would in- 
crease by a fraction {SVb)/V where SVb is the increase 
in volume occupied by the bubble and V is the volume 
occupied initially by Phase 1. This can be made small by 
choosing a large container. We can also imagine a situa- 
tion where the transition from Phase 1 to Phase 2, does 
not involve a significant change in density. This might 
happen, for example, in the case of a fluid to superfluid 
transition. The bubble wall then just marks the location 
where Phase 1 is being converted to Phase 2, without the 
wall actually exerting any forces on the material outside. 
In this case, there is no transport of material involved. 



In fact, the two phases need not even be fluid, they can 
be solid, say the normal and superconducting phases of a 
Type I superconductor. Our only requirement is that we 
have two different sound velocities and a phase bound- 
ary velocity that satisfies ci > v > C2- If the velocities of 
ordinary (compressional) sound do not satisfy these con- 
ditions, one could consider some other excitations. How- 
ever, it is advantageous to consider the excitations with 
the largest velocities that satisfy the conditions since then 
the Hawking temperature is highest. 

In this setup, the Hawking temperature can be rela- 
tively large because the gradient of Cs — v occurs within 
the phase boundary whose thickness will be of micro- 
scopic dimensions, given by the inter-molecular spacing, 
or the coherence length of Cooper pairs, or some other 
such length scale. The sound velocity depends on what 
kind of excitations we are considering. We would like 
to take a large sound velocity, since this gives a larger 
Hawking temperature. With ordinary sound (compres- 
sional waves) typical sound speeds are ^ 300 m/s. This 
leads to the estimate: 

where Scg is the change in the sound velocity across the 
phase boundary and ^ denotes the thickness of the phase 
boundary. This is quite large compared to Hawking tem- 
perature estimates in earlier proposals. The thermal fre- 
quency is V = kBT/h ~ 1 GHz and the power emitted is 
(TsT^fj ^ 10^ pW/cm^ where the sonic Stefan-Boltzmann 
constant CTs = 7r^fc|;/60cf involves the sound speed and 
two polarizations have been assumed |^. 

Now let us consider some specific systems and varia- 
tions on the general scheme described above. The first 
point to note is that we must necessarily work at low 
temperature since the Hawking temperature is of order 
0.04 K. All substances except the isotopes of Helium, so- 
lidify at such a temperature. So we need either to work 
with a solid, or with Helium. If we choose ''He, we will 
be working with a superfluid. '^He however need not be 
in the superfluid phase. 

An existence proof that the conditions Ci > v > C2 can 
be met can be found by considering the Abelian Higgs 
model in a cosmological setting. (The non-relativistic 
version of the Abelian Higgs model is the Ginzburg- 
Landau model used to describe simple superconductors.) 
For a certain set of parameters, the Abelian Higgs model 
is known to undergo a first order phase transition, cor- 
responding to Type I superconductors. Imagine that the 
system is supercooled, so that it is stuck in the metastable 
state labeled Phase 1 in Fig. |2 Excitations in this 
metastable state are massless gauge particles (photons) 
and massive scalars. If these excitations are in equilib- 
rium, they form a relativistic fluid for which the equation 
of state is p = p/3 (in units where the speed of light is 
one). So the sound velocity in the supercooled state is 
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FIG. 2: A sketch of the free energy V of the order parameter 
(j). The force driving the bubble expansion comes from the 
difference in energy e between the two local minima. The 
drag on the bubble arises due to scattering of the excitations 
surrounding the bubble with the bubble wall. 



ci ~ 1/ \/3. If there was no plasma outside the bub- 
ble, the bubble wall velocity would be the speed of light. 
However, the plasma exerts a drag on the bubble wall, 
and reduces its speed somewhat, and one has v ci. (In 
the similar case of the cosmological electroweak phase 
transition, calculations yield v 0. 1 Hi.) Inside the 
bubble, in Phase 2, all the excitations are massive and 
non-relativistic. Now the equation of state is p ~ and 
the speed of sound C2 is non-relativistic, say 300 m/s 
which is much less than v. Therefore the conditions 
ci > V > C2 are met, and the bubble should emit sonic 
Hawking radiation. Also note that any heat released by 
the phase transition within the bubble does not propa- 
gate to the outside since the bubble wall is moving at 
supersonic speeds. So oscillations of the order parameter 
in the true vacuum are inconsequential. The only danger 
arises from excitations living outside the bubble {i.e. or- 
der parameter oscillations in the metastable state) that 
might scatter off the bubble wall. However, if the system 
is sufficiently supercooled, the density of excitations in 
the metastable state is expected to be very low. 

A one dimensional version of our setup consists of a 
cylindrical container filled with Phase 1 (Fig. O . A per- 
turbation applied on the right-hand end of the cylinder 
will trigger bubble nucleation of Phase 2. The bubble of 
Phase 2 will start growing and the phase boundary (or 
bubble wall) will propagate to the left. Unlike the bubble 
of Fig. ^ the area of the sonic horizon does not change 
with time and hence the setup is easier to analyze. The 
disadvantage of the one dimensional setup is that there 
is more interaction of the system with the walls of the 
container and hence there is greater danger of instabili- 
ties. In particular, heat released on the Phase 2 side can 
be transmitted to the Phase 1 side through the walls of 
the container. 

A one dimensional condensed matter situation where a 
phase boundary propagates at supersonic speeds occurs 
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FIG. 3: A one dimensional version of the experimental setup 
considered in this paper. The system is kept at constant tem- 
perature and pressure. 



in '^He. The boundary of a bubble of A phase inside a 
metastable B phase has been experimentally seen Ti^ to 
propagate at speeds of up to 67 cm/s, while fermionic 
quasiparticles excitations in '^He-A propagate at only 3 
cm/s. On the ^He-B side, these excitations are can prop- 
agate at 55 m/s and so the setting seems right for the 
creation of a sonic horizon. However, note that the bub- 
ble wall velocity is in the wrong direction in this case 
- we would like the wall of an A phase bubble, where 
the sound velocity is small, to propagate at high velocity 
into the B phase where the sound velocity is high. It re- 
mains to be seen if the A-B phase boundary can provide 
measurable Hawking radiation Q (also see Sec. 29.3 of 
111). 

A novel feature of our proposal is that it is an example 
of a "non-driven" sonic horizon. In other examples, one 
has to arrange external forces to produce a sonic horizon. 
For example, in fluid flow, the fluid has to be driven un- 
der pressure. If the external conditions are designed to 
maintain a sonic horizon, the Hawking radiation backre- 
action on the sonic horizon will be lost. In the bubble 
nucleation proposal, however, the procedure is to nucle- 
ate a bubble of Phase 2 and then leave the system to 
its own devices. If the propagating phase boundary does 
emit Hawking radiation, it will experience a backreac- 
tion that will slow it down. Eventually the bubble wall 
velocity will become subsonic with respect to the sound 
velocity in Phase 2. In other words, the prediction is that 
the condition ci > v > C2 is quantum mechanically un- 
stable and eventually the system will go to a state with 
ci > C2 > V. This is no different from the evaporation of 
black holes due to Hawking radiation. 

To estimate the evaporation time, note that the power 
in the emitted Hawking radiation is: ^ asT^^^A where 
A ^ AttR'^ is the area of the horizon when the bubble 
radius is i? = vt. The energy lost as a function of time 
is therefore: Eiost ~ cr sTf^^ivv'^t^ . The kinetic energy 
of the phase boundary in the absence of the losses would 



be: Ekin ~ SAd^ where S denotes the surface energy 
density. Therefore the evaporation time scale is given by 
equating the energy lost to the total kinetic energy: 
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where / = v/ci. The evaporation time depends very 
crucially on the surface energy density. We would like 
to have a large surface density so that the evaporation 
is relatively slow. The surface density depends on the 
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system that is being used. However, we can get an idea 
for its value by considering the phase boundary to be a 
domain wall. For a domain wall in a X4>'^ model, we have: 

Hence if the coupling constant A is small, the evaporation 
time will be long. This is exactly as for black holes where 
the Hawking evaporation time is long because gravity is 
so weakly coupled. One way in which black hole evapo- 
ration is different from dumbhole evaporation is that the 
black hole necessarily gets hotter as it gets smaller. The 
dumbhole may get hotter or colder with evaporation, de- 
pending on the gradient of Cs{t,x.) at the location of the 
sonic horizon v = Cs(t, x). 

If it is difhcult to find a system in which the condi- 
tion V > C2 holds, one could consider driving the phase 
boundary by some external forces. This might be easy 
to implement if the substance undergoes a phase transi- 
tion when an electromagnetic field is applied. Then the 
electromagnetic fields could be varied externally to drive 
the phase boundary at high velocities. The only danger 
is that the application of dynamical external forces could 
give rise to additional heating that could swamp the sig- 
nal. We will not consider the "driven" setup any further 
here. 

We now discuss potential difficulties in implementing 
the current proposal. 

The first danger is that the phase boundary may have 
instabilities during propagation. Such instabilities need 
to be avoided since they will introduce noise in the 
Hawking emission. Fingering and other instabilities are 
system-dependent, and the success of the proposal will 
depend on finding a system in which such instabilities 
are absent. 

An important issue is that, since the phase boundary is 
propagating, eventually the whole container will be filled 
with the second phase. So there is only a limited time 
{L/v — L/ci) ^ L/v, where L is the size of the container, 
to make measurements. This time had better be long 
enough to detect excitations at the thermal frequency 
v — ksTsH/h ^ 1 GHz and we obtain a lower bound on 
the required length of the container: L > hv / {ksTsu) ^ 
10~^ cm. This condition does not seem difficult to meet 
in the laboratory where samples are typically on the scale 
of centimeters. However, a sample of 1 cm size will only 
provide ~ 10~^ s to make measurements. 

Another potential difficulty is that Phase 1 is 
metastable and might decay spontaneously. Hence it has 
to be protected from perturbations. However, the Hawk- 
ing radiation itself will perturb Phase 1 and the concern 
is if the system might "self-destruct" . To avoid this pos- 
sibility, Phase 1 should be immune to sonic perturbations 
at the Hawking temperature. 

The bubble will have perturbations that are produced 
during the nucleation process. Even quantum nucleation 



leads to fluctuations on the bubble |l5l llq . The ampli- 
tude of fluctuations is suppressed if the bubble surface 
density, S, is large. Hence the effects of bubble pertur- 
bations can be reduced, as well as the evaporation time 
made large, by choosing a system with a large S. 

In summary, we have shown that propagating phase 
boundaries can be sonic horizons provided the propa- 
gation velocity satisfies the condition ci > v > C2. In 
systems where this condition is met, the Hawking tem- 
perature can be relatively large, leading to the exciting 
possibility that experiments might detect Hawking radi- 
ation. In addition, in these systems the sonic Hawking 
radiation will lead to evaporation of the sonic horizon, 
providing a yet closer analog of gravitational black holes. 

I am grateful to Ana Achucarro, Dan Akerib, Craig 
Copi, Arnie Dahm, Harsh Mathur, Rolfe Petschek, Bill 
Unruh, and Grisha Volovik for discussions and com- 
ments. This work was supported by DOE grant number 
DEFG0295ER40898 at Case. 



[1] W.G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 

[2] S.W. Hawking, Nature 248, 30 (1974); Comm. Math. 

Phys. 43, 199 (1975). 
[3] W.G. Unruh, in "Artificial black holes", eds. M. Novello, 

M. Visser, and G. Volovik, World Scientific, River Edge, 

USA (2002). 

[4] L.J. Garay, J.R. Anglin, J.I. Cirac and P. ZoUer, Phys. 

Rev. Lett. 85, 4643 (2000). 
[5] L.J. Garay, J.R. Anglin, J.I. Cirac and P. ZoUer, Phys. 

Rev. A63, 023611 (2001). 
[6] P.O. Fedichev and U.R. Fischer, Phys. Rev. Lett. 91, 

240407 (2003). 
[7] M. Visser, Class. Quant. Grav. 15, 1767 (1998). 
[8] T.A. Jacobson and G.E. Volovik, Phys. Rev. D58, 

064021 (1998). 

[9] W.G. Unruh and R. Schutzhold Phys. Rev. D68, 024008 
(2003). 

[10] The sound velocity that varies within the domain wall 
is not for the usual compressional waves but for certain 
excitations that are very specific to the '^He-A order pa- 
rameter. 

[11] To connect with the Unruh setup, it is more convenient to 
transform to the rest frame of the phase boundary. In this 
frame the fluid is fiowing with velocity v which is subsonic 
outside the bubble and supersonic inside. Note that the 
process of cavitation, in which the bubble wall simply 
pushes the fluid and leaves behind an empty cavity, does 
not provide a sonic horizon. In this process, in the rest 
frame of the bubble wall, the fluid is also at rest. 

[12] G.D. Moore, JHEP 0003, 006 (2000). 

[13] G.E. Volovik, "The Universe in a Helium droplet", Ox- 
ford University Press (2003). 

[14] D.S. Buchanan, G.W. Swift and J.C. Wheatley, Phys. 
Rev. Lett. 57, 341 (1986). 

[15] T. Vachaspati and A. Vilenkin, Phys. Rev. D43, 3846 
(1991). 

[16] J. Garriga and A. Vilenkin, Phys. Rev. D45, 3469 (1992). 



